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^ ■ Abstract 
O . 

. Starting with the Dirac equation outside the event horizon of a non-extreme Kerr black hole, we 

develop a time-dependent scattering theory for massive Dirac particles. The explicit computation 
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of the modified wave operators at infinity is done by implementing a time-dependent logarithmic 
phase shift from the free dynamics to offset the long range term in the full Hamiltonian due to the 
presence of the gravitational force. Analytical expressions for the wave operators are also given. 
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I. INTRODUCTION 



In this paper we develop a time-dependent scattering theory for massive Dirac particles 
outside the event horizon of a non-extreme Kerr black hole manifold. Such a solution of the 
vacuum Einstein field equations discovered by Kerr in 1963 describes an asymptotically flat 
space-time containing nothing but an eternal, rotating black hole, and has been generalized 
to the case of charged, spinning black holes by Newman et al. in 1965. Although it is not 
the most general model of the exterior region of a black hole we can analyze theoretically, 
it represents indeed the most realistic model in astrophysics since in general black holes are 
embedded in environments that are rich in gas and plasma and, consequently any net charge 
is neutrahzed by the ambient plasma (see for instance Misner et al. (1999) p. 885). 
The aim of scattering theory on curved space-times is to provide a detailed description of 
the asymptotic behavior in time of some field (in our case Dirac fields). A general feature 
of scattering theories on black hole manifolds is that two asymptotic regions are present. In 
fact, in the analysis of the propagation of fields outside the event horizon of a black hole we 
commonly adopt the point of view of an observer static at space-like infinity, who in turns 
perceives the event horizon as an asymptotic region. 

There are mainly three motivations for the analysis of our problem. The first one, quite 
general on its own, aims to provide a deeper insight into the physics of black holes with 
particular attention to the quantum field theory on curved space-times. The second one 
is the rigorous mathematical analysis of quantum effects in general relativity such as the 
Hawking effect (see Hawking (1975), Bachelot (1997), (1999), (2000), and Melnyk (2004)). 
Finally, the third motivation is dictated by the study of resonances, i.e. of the complex 
frequencies which are poles of the analytic continuation of the scattering operator. For such 
studies in the Schwarzschild geometry we refer to Bachelot, and Motet-Bachelot (1993), and 
Sa Barreto, and Zworski (1998). 

Concerning the time-dependent scattering theory for Dirac particles in a Coulomb field we 
find some early publications by DoUard (1964), by DoUard, and Velo (1966), and by Enss, 
and Thaller (1986). The time-dependent scattering theory for classical, and quantum scalar 
fields on the Schwarzschild metric was first obtained by Dimock (1985), and by Dimock, and 
Kay (1986a,b), (1987). In the same geometry Bachelot developed the scattering theory for 
electromagnetic fields (1991), and Klein-Gordon fields (1994). Regarding Dirac fields in the 
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Schwarzschild geometry, Nicolas (1995a) presented a scattering theory for classical massless 
Dirac particles, and Jin (1998) constructed wave operators, classical at the event horizon and 
DoUard-modified at infinity, in the massive case. Moreover, Melnyk (2003) gave a complete 
scattering theory for massive charged Dirac fields in the Reissner-Nordstr0m metric. Finally, 
Daude (2004) proved the existence and asymptotic completeness of wave operators, classical 
at the event horizon, and DoUard-modified at infinity, for classical massive Dirac particles 
in the Kerr-Newman geometry by means of the Mourre theory (see Mourre (1981)). For the 
nonlinear Klein-Gordon equation on Schwarzschild like metrics partial scattering results by 
means of conformal methods have been obtained by Nicolas (1995b). A complete scattering 
theory for the wave equation, on stationary, asymptotically flat space-times, was developed 
by Hafner (2001). 

Whenever we attempt to analyze the scattering properties of fields outside the event horizon 
of a Kerr black hole, we are faced with some difficulties which are not present in the picture 
of the Schwarzschild metric. First of all, the Kerr solution is only axially symmetric since it 
possesses only two commuting Killing vector fields, namely the time coordinate vector field 
dt and the longitude coordinate vector field d^p. This implies that there is no decomposition 
in spin-weighted spherical harmonics. Moreover, another apparent difficulty is due to fact 
that it is impossible to find a Killing vector field which is time-like everywhere outside the 
black hole. In fact dt becomes space-like in the ergo-sphere, a toroidal region around the 
horizon. This implies that for field equations describing particles of integer spin (wave equa- 
tion, Klein-Gordon, Maxwell) there exists no positive definite conserved energy. For field 
equations describing particles with half-integer spin (Weyl, Dirac) we can find a conserved 
L2 norm with the usual interpretation of a conserved charge. Hence, the absence of station- 
arity in the Kerr metric is not really a difficulty for the scattering theory of classical Dirac 
fields. Nevertheless there are only few analytical studies of the propagation of fields outside 
Kerr black holes. 

Our work represents a new approach to the results obtained by Jin (1998), and by Hafner, 
and Nicolas (2004). Firstly, our method is based on an integral representation for the Dirac 
propagator outside the event horizon of a non-extreme Kerr manifold: this is new in this 
context. Secondly, we are able to compute explicitly the wave operators (Dollard-modified) 
at infinity. Moreover, by computing the wave operators (classical) at the event horizon and 
introducing suitable global wave operators, it should be possible to give an alternative proof 
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for the asymptotic completeness to that one presented by Hafner, and Nicolas, and to cal- 
culate the scattering matrix. This will be done in the next future. 

Let us briefly describe the contents of this paper. In Section IH] we give the integral rep- 
resentation for the Dirac propagator in the exterior region of a Kerr manifold, and the 
asymptotic behavior for the radial solutions. These results are essential to the development 
of our theory. In Section UTTl we define the free dynamics asymptotically at infinity, and we 
introduce the so-called DoUard- modified wave operators. In Section IIVI we compute explic- 
itly the phase shift we need to implement in the free dynamics. In this section the main 
result is Theorem IIV.II where we give an integral representation for the DoUard-modified 
wave operators. 



II. PRELIMINARIES 



In Boyer-Lindquist coordinates {t,r,'d,ip) with r>0, 0<'i9<7r, 0<(y9<2'7r the Kerr 
metric is given by (e.g. Wald) 

, n A - a^sin^-i? , n 2asin^'(9(r^ + - A) , , E , ^ ^ , , , ^.^E , 9 

ds^ = df + ^ -dtdip - —dr^ - T.d§^ - (r^ + a^f-d^^ (1) 

E E A E 

with 

E := E(r, 9) = r'^ + c? cos^ A := A(r) = - 2Mr + 

and 

E := E(r, i?) = 1 — a^7^(r) sin^ 7(r) 



where M, and a are the mass, and the angular momentum per unit mass of the black hole, 
respectively. Here, a is allowed to be zero, so that our results apply also to the Schwarzschild 
metric. Moreover, we will always work in the non-extreme case > which implies that 
the function A has two distinct zeros at the Cauchy horizon tq = M — \JM'^ — a^, and at 



the event horizon ri = M + \/ M'^ — . Notice that A > for r > r^. Outside the event 
horizon of a non-extreme Kerr manifold the integral representation of the Dirac propagator 
for a particle of mass me, charge e, and energy lo is (Theorem 5.4 in Batic, and Schmid 
(2006)) 

/+00 
du^e-' E x = [u,{},^) (2) 



4 



where (-j-) is the positive definite scalar product 



+00 



du / d(cost9) / d(pV'£ i/j^ -00, 



(3) 



and -00 denotes some initial data in C^(0)^ with Q, :— x S'^. Here, m e E is the so- 
called tortoise coordinate defined by du/dr — (r^ + a^)/A. Notice that the variable u has 
the property to approach —00 as r approaches the event horizon. Furthermore, A; e Z is 
the azimuthal quantum number, j e Z\{0} labels the eigenvalues of the angular operator 
arising from the Chandrasekhar separation of the Dirac equation into a radial, and an 
angular system of ODEs (see Chandrasekhar (1976)), and according to the Chandrasekhar 
ansatz -0^-^ has the following form 



R. 



kj 



U 



{RtU^)s'J,Mj 



,i{k+^)<P 



(4) 



with Rj' = {bI^_^R^^_^Y'^ and S^^ = {S^-'_, Sly^_^_)^ the radial, and angular components of 
the spinor ipj', respectively. For further properties on the angular eigenfunctions Sj^ we 
refer to Finster et al. (2000), and to Batic et al. (2005). As shown in Lemma 6.1-2 in Batic 
and, Schmid (2006) in the limit u — > —00 the radial solutions behave for w < i^i < as 
follows 



R'Jiu) 



(5) 



with 



^(o)g^iQoui guch that |cL°^|2 



+ c. 



(0)|2 

■+ I 



7^ 0, and 



< d — 4ac+, k+ — 



ri - To 



rt + a^ ' ■ ' ■ 2{rf + a^) 

where is the surface gravity at r = ri. For u +00 wc have to distinguish between 
the cases |a;| > nie, and |a;| < rUe- In the first case the asymptotic behavior of the radial 



functions is ior u> uq > 



< 



(6) 
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where = c±eT^*("o) with c± e C such that |c_P + |c+| V 0, 

coshO sinhO \ 1 /c^ + me\ , Mm^ 

sinhe coshe / 4 v^-^e/ « 



K = e{u)^/u!'^ — ml, e(a;) — 



+ 1 if > +m, 



e 



-1 if a; < —rrie 



and 



with 



K 2M{\j{uj)-2mia') MmlPkj{uj) M^m 



e 



4 \ ^4 ^6 

Notice that for some e > R>{u) is smooth in a; e (— 00, — mg — e] U [me + e, +00). 
For |tc;| < rrie there are two fundamental solutions with exponential decay, and growth, 
respectively which are smooth in a; G [—me + e, me — e] . In order to disregard the unphysical 
solution with exponential growth we normalize the radial solutions i?< (u) by imposing that 
R^{u) ^ 1 as — > 00. 



III. DOLLARD-MODIFIED WAVE OPERATORS 

Because of the presence of two asymptotic regions ^ ±00) we need to specify for each 
of them an asymptotic dynamics. For u +00 we define the free dynamics by replacing the 
Hamiltonian (4.2) in Batic, and Schmid (2006) by its formal limit H^o when M — > 0. Notice 
that in this case the Kerr metric goes over to the Minkowski metric in oblate spheroidal 
coordinates (OSC). Since du/dr = 1 for M = 0, and r can be extended to negative values, 
we can identify the tortoise coordinate u with the spatial variable r. We consider Hoo as an 
operator acting on the Hilbert space 



"Hoo = L2 (fi, rf/Xoo)^ , djioo = \l , ^ ^ du d(cosi?) dip. 
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Moreover, the Hamiltonian Hoo is formally self-adjoint with respect to the positive scalar 
product (Sec. 7 ibid.) 

Jo. 

it is essentially self-adjoint on C^(f2, rf/ioo)'^, and it has a unique self-adjoint extension on 
the Sobolev space W^'''^{VL^d^^Y . The Dirac propagator for a particle of mass me, charge 
e, and energy uj in the Minkowski metric expressed in OSC is (Lemma 7.2 ibid.) 



where ^ Co^(fi, d/ioo)', a(iJo 



(— oo, — rrie] U [me, +cxd), and 

/i?^^:^°°)(«)5^f_(^)\ 



R, 



(9) 



Notice that the angular components of the spinor ipt'''^'^^ are the same as those for tpjy since 
the angular operator arising from the Chandrasekhar separation of the Dirac equation into 
a radial, and an angular system of ODEs does not contain the mass parameter M. The 
asymptotic behavior of the radial functions is (Lemma 7.3 ibid.) 



Tiuo) 



u)] 



kj,{oo) 

±,UJ 



{u)\< 



(10) 

where /^°°^ = c±e^'"^'^° with c_, c+ G C such that |c_p + |c+p ^ 0. Moreover, 9, and k are 



defined as in Section HTl and 



Pkj{uj] 



2k ^''"-^ ■ ■ ^K^ 

Clearly, Rj''^°°\u) is smooth in G (— cxo, — mg — Jl] U [rrie + Jl, +oo) with Jl > 0. In 
what follows we consider the full Hamiltonian H as an operator acting on the Hilbert space 
Ti. = L2{^,dfi)'^ with dfi = \/t, du d{cos§) dip. Since we want to compare H with Hoo we 
define the smooth bounded identifying operator Too : ?ioo — ^ 'H as follows 



(11) 



with cut-off function Xoo ^ C°°(R) such that 



Xoo{u) 



1 if u > Mo 

(12) 

if u < Mo 



with < r]oo{u) < 1 for Mq > u > Mq > 0. Because of the long-range nature of the 
gravitational potential we define asymptotically at infinity the following Dollard-modified 
wave operators (see DoUard (1964)) 

Wt^.^t^ = s- hm e-^*Tooe^^-*e^^W4°^\ 4°°^ G Ho. (13) 



t— >±oo 



under the conditions that the phase shift operator 6{t) commutes with H^, and that the 
limit exists in the strong sense in TC. If we apply Q, and (jH)) we can express (fT^ as follows 

E E^c.'''(^)V^^^'^°°Ha^)(^^^°^Vr^)(oo). (14) 
jez\{o} fcGZ 

In preparation of the next results we define the intervals B^{±me) := {to' G ]R| |u;'±me| < e}, 
/+ := {uj G (7(-ffoo)|0 < u — rrie < /!}, and /_ := {uj G cr(i/oo)| — /U < + mg < 0} for some 
e > 0. The following lemma controls uniformly in t the contributions of the frequency sets 
B^l^rrie), and J± to the wave operators asymptotically at infinity. 

Lemma III.l Let E : M — > B{H), and E°° : a{Hoo) — ^ B{Hoo) be the spectral families 
associated to the self-adjoint operators H and Hoo, respectively. In particular, let us consider 
their restrictions on the Borel sets := B^{—mf.) U B^inif,), and / := /+ U i.e. : 
B, — > Bin), E^ -.1 — > BiUoo) with 



E.:= j XBX^')dE{J), Ef:= j 



Xi{uj) dE^ico) 

where Xb^, o-nd xi o'^e the characteristic functions of the intervals B^, and I. Then, for 
every k > 0, and tp^""^ G C^{^, dn^oY there exist constants Jl, e > such that for every t 

\\Uoo{t)4'^^ - e-^^*(/rf - E,)Jooe^^-*(/rfoo - i?2°)e^^W4°°^lk < ^ 
where Uoo{t) ■= e~^^^Iooe^^°°^e''^^^\ and 5{t) is a phase shift operator commuting with H^o. 



Proof. Since there is no risk of confusion we omit to write explicitly the superscript (oo) 
attached to the initial data. Let us define 



By adding and subtracting the term e~*^*Xooe*''^°°*(/(ioo — 6**^*^*^0 to AW we obtain 
AW = e-'^^Iooc'^^^e'^^^^E^iPo + e~'"*E,I^e'"^*e'^^^^iPo - e-'^^E^Io^e'^'^^E^e'^^^^iPo- 

Since e"'-'^*, and e*^°°* are unitary, I^o is bounded, e*''*^*) commutes with E-', and E^ com- 
mutes with e~*^*, we have 

\\AW\\n < 2\\E^^|;o\\n^ + \\E,U^{t)Mn- (15) 

Let us recall that 



Without loss of generality we consider the above expression for the set /+. Taking into 
account that the idempotent projector is hermitian, and making use of (5.8) in Batic, 

and Schmid (2006), it follows that 

= J du;xi+{oj)-^{ipo\E°°{uj)ij;o) = J dujxi+{uj){'ipo,u,\'tpo,u,)i,oo 

where V^o,a; = (^^o)(^) ^ f)oo is the representative of the element ^ ^oo, the map 
JF : Hoo — > ^oo is unitary, and ('?/'o,u)|'«/'o,a))[)oo = ["^o.ujp- For the definition of i)^ see 
Theorem 5.2 (ibid.). Thus, we obtain 

\\E^Mn^= j dujxi{uj)\i>0,.? ■ 

Since the r.h.s. of the above expression converges to zero with the length of /, by choosing 
/I sufficiently small we can arrange that 

WE^Uh^ < '^/4. (16) 

According to Daude (2004) the strong hmit of Uoo exists. Moreover, Lemma 5.3 in Batic and 
Schmid (2006) imphes that the point spectrum of H is empty. Hence, for a fixed e it 
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results that for e 0. Such a hmit is uniform for G C C 7i with C compact. 

Since Uooit)ipo is continuous in t and converges for t ±oo, the closure of {t/oo(^)'?/'o} is a 
compact set. Thus, for a given A; > we can always choose e sufficiently small such that for 
every t G M 

WE^U^moWn < \- (17) 
The estimates (ITHjl . and (fTTj) together with (fTHjl complete the proof. □ 

Let us now introduce the set A = {uj' E W\ — ojq < uj' < oJq] for some ujq > rrie. In the 
following remark we explain how to bring the limit t ±oo inside the integral over u' , and 
why we can consider a finite number of terms in the sums over j', k\ j, k entering in (jl4p . 

Remark III. 2 Let E be the spectral family defined in Lemma \III. 11 and let us consider its 
restriction on the set A 

: A ^ B{n), := j xa(^') dE{J) 

where xk is the characteristic function of the interval A. Since Uoo(t) converges strongly to 
fort -> ±oo (see Daude (2004)), and ^Af/oo(t) converges to E^Wf'^^ for t ±oo it 
follows that 

Finally, taking into account that the spectrum of H is purely absolutely continuous, and that 
Ex converges strongly to the identity for ujq — > +oo it results that E^^W^^ converges strongly 
to for ujq +00. By means of a similar argument it can be shown that we can restrict 

our attention to a finite number of quantum numbers j' , k' , j, k. 

IV. THE WAVE OPERATORS ASYMPTOTICALLY AT INFINITY 

Let us define the sets Qj := [—uJo, —nie — e] U [nie + e, Uq], Qji := [— mg + e, mg — e], and 
Qui := (— cxo, —rUe — Ji\U [nie + Ji, +oo). Taking into account that admits the following 
asymptotic expansion in the spatial variable u 

VE= l + 7^(u,^9), |7^(M,^9)|< — 

if we apply Lemma IlII.ll Remark IIII.2t and if we make use of (jH), pOj) . and Lemma 
6.1 in Batic, and Schmid (2006), concerning the asymptotic behavior of the radial functions 
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R^{u), then we end up after a lengthy computation with the following expression for the 



DoUard-modified wave operator (jT4| 



with 



'^'^(oo),//7/V^0 



duj' Yl ^5 jiToo / "^""^ 



'■I/II 



\i'\<j'o I-jI--?" I'^I^'^o 



du e'("-"')*+*^(*)Fj^//,(cu', cu, m), (18) 



fijij 



UJ 



2t: 



+e 



+ h'^-' {uj , UJ, u) 

/^H, (19) 



<^(u) \ —inu ■ 



\ f- f- + 



U 



+ T^" {U ,UJ,U) 



+ 



+e^^"A'(u;',cj) 



^kn' \f-fr + 



U 



+ Tj^^ {UO ,00, U) 



g'^iu) (20) 



where g''^{uj) = {ip!y'^°°'^\^Po°°^)(oo), and 



Z{uj , uj) 
X{uJ , uj) 
Ziio , lS) 
X{lo , iS) 

Furthermore, 



cosh 6(c<j ) cosh Q{u!) + sinh Q{u ) sinh Q^uj), 

cosh 9(^7 ) sinh 9(u;) + sinh 9(ci; ) cosh Q{uj), 

sinh 9 (tu) + r(co' ) cosh 9(^7), t{uj ) = — 

rrie 

cosh 9(^7) + t{uj ) sinh 9(ci;), 



UJ + i\ mi — UJ 



rf(cos^) [s'j^_{^)Stl.m + ^5;+(^)^^f+(^) 
rf(cos^) sin^^ fe'_(^)^5-(^) +^1' +(^)^'?+(^) ) ' 
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and {u',uj,u), h!^-' {uo\uj,u) are 0{u^'^) with 




\h^-' [u ,u,u)\< , \h^-' [uj ,uj,u)\< 



where 



^kj j 



n(2).(°o) 



Dfir\^)+Tf^],.(J)^D%.(.)D 



T,kj 



Notice that 2t_|_"'"' , ^B^"''' , h'^-' , /i^-'-' G C°°(n7 x Qm). Similar formulae and considera- 
tions hold also for 21^-'"' , , r^-'-' , and r^-'"' . Moreover, since the angular eigenfunctions 
S^, _^{^), S^-'j.{'d) are smooth in u' , and u (Sec. II in Batic et al. (2005)), it follows that 
Aj^jj', Bf^--i G C^iVLj X VLiii)- Finally, by means of the Holder inequality, and taking into 
account that the angular eigenfunctions are normalized it can be easily shown that ^kjj'i 
B^jj' < 2. Although (HH)), and look quite complicated, we shall show in this section 
how to reduce ()18|) to a more amenable form. 

To determine which phase shift operator we have to introduce in order that ()13|) makes sense 
we can proceed as follows. Let us consider for a moment the classical definition of the wave 
operators asymptotically at infinity. If we proceed exactly as in the present section, we end 
up with an expression analogous to (fTS|) but without phase shift whereas (fTn|) . and ()2()|1 
remain unchanged. At this point we just need to analyze for u G [mq, +oo) the long-time 
behavior of the integrals containing terms of zeroth order in the expansion in powers of 
Lemma fA.2l suggests the following definition for the phase shift, namely 



5[t) ■= -a{uj) Log i-t 



Mml 

a{ijj) := e{uj)- 



(21) 



where as in Dimock, and Kay (1986a) Log t is defined for alH 7^ by Log t := sgn(t) log \t\ 
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Theorem IV. 1 Let W^-^ be defined as in (fT^ with phase shift 5{t) specified by (j^H). Then 
for every ?/'^°°-' G C^{Q,dfioo)'^, 




Proof. Let us begin with the foUowing observations. The function g^^{uj) = 
{ipj''^'^^\ip()^'^){oo) entering in (fT^ . and (pUj) belongs to the space of smooth rapidly de- 
creasing functions since by going from oblate spheroidal coordinates to Cartesian 
coordinates g'^^ is simply the Fourier transform of smooth initial data with compact support. 
Moreover, for to' E ^i/ii by means of (7), and (8) in Batic et al. (2005) it can be checked 
that Qi'^'' , QS^-'"' , h'^-' , and h'^-' together with their first cj-derivatives are polynomially 
bounded in u which implies that their products with g^^ are again Schwarzian functions. 
Analogous considerations hold for Sl^-'-' , QS^^"* , r^"'-' , and t^^-' . Finally, notice that 

Z{u! ,uj ) = — , X{uj , —UJ ) = 0. 

Theorem IB. 81 ensures that for uj' G fi// there is no contribution to the wave operators 
asymptotically at infinity. Concerning the case uj' G fi/ we apply Theorems IB.3| IB.4| IB.5[ 
and IB. 61 to obtain 

/, p+oo , 

E E / duu^- 

j duj Z{J,uj)A^.^>{J,u)g^\u)J^f^e'^^-'''^'+^'^'^-^>+'^^^^ 
The result follows by applying Theorem IB. 31 □ 
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APPENDIX A: OSCILLATING INTEGRALS 

Let LU G fi///, to' G fi/, K, be defined as in Section |Tll k = e{uo' ) \J uo'"^ — ml, uq > 0, 
and a' := Mml/n . Moreover, let JF g where ^ denotes the Schwarzian space of 

smooth rapidly decreasing functions. 

Lemma A.l 



Ji± := lim / duu^"^ \ dcjj-(cj)e'("-" = ij J^(W) = 0, (Al) 

Qui 

r+°° if II 

J2± := lim / duu^''' / rfcjj^(cj)e^('^-'^ +'^)" = if J^{-J)={] (A2) 



' "0 

Qui 



where the subscript ± attached to I1/2 corresponds to the ± entering in the exponents of the 
integrands. 

Proof. Concerning Xi j-, after the introduction of a convergence generating factor e"*^" with 
0" > we can apply Fubini theorem to obtain 

/, r+co , 
du J^{uj) e'^^-^ / dr r±^° e""*" (A3) 
Ji 
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with r = m/mo, and 

i± := 'z±{uj) = [cr =1= i{K — k)]uq. (A4) 

By means of 13.2.6, and 13.1.3 in Abramowitz, and Stegun we can perform the integration 
over r, and Xi ± becomes 



^ = ul^'^Vil ± la) lim lim I duo T{uj)z e^^""'^')* 

Qui 

y lim lim [ diuJ^(uj)e^'^M(l,2±ia';z±)e'^'^-'''^' (A5) 

l±ia t-*±oocr^o+ J 

Qui 

where M denotes the confluent hypergeometric function of the first kind. Regarding the 
second integral in the above expression, let us define the function 

J^^[uj) := J^(cu)e-^±M(1, 2±ia'; i±)e^('"-"')*. 
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Using the integral representation for the confluent hypergeometric function of the first kind 
13.2.1 (ibid.), we obtain the estimate |M(1, 2±ia'] < \/ 1 + a'^ e'^"° from which it follows 
that |jFo.(co')| < \/ 1 + a''^\J-'{uj)\. Hence, the Lebesgue dominated convergence theorem can 
be applied to the second term in ()A5|) which simplifies to 



lim / du (u) M {1,2 ±ia']z±) e 
-^±00 J 



i(uj—u} )t 



Qui 



with J^(co') = jF(co')e ^±,and£j- := ^iui{K' —k). Since F(ci;) G and M{l,2±ia' ;'z±) 

is bounded by + a , the Riemann-Lebesgue lemma implies that the above integral is 
zero for t ±00, and (IA5I) becomes 



Ji± = r(l±m) lim lim f du ,e^(^-^')*. (A6) 

Let us rewrite the fraction entering in the above integral as follows 



Notice that \{k — t^)/[(y ± — < 1. Moreover, in the case t{uo') = — e(u;) the function 
T{uj)/ (k — k) is integrable, and we can immediately apply the Lebesgue dominated conver- 
gence theorem to take the limit a ^ 0^ inside the integral in ()A6|1 whereas for e{u') = e(u) 
we make in the above expression the substitution 



K — K UO — UO' K — K '' 



and observe that T{uj)/{uo — J) has rapid decay for \uo\ — > 00, and it is integrable since 
J^{oj') = 0. Hence, ()A6|) reduces to 

Finally, the Riemann-Lebesgue lemma implies that Ii^± = 0. We omit the proof for 22,± 
since it resembles that one for Ti±. □ 
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Lemma A. 2 LetXi^±, andl2^± he defined as in Lemma \A.l\ Then 



if t ^ —00 



1 +00 

'^T{J) e-*"'i°g(-*) ift^-oo ' 

z/ 1 ^ +00 

^J^(-Cj') 6-2^"^'* if -00 ' 

if t ^ —00 



Proof. We compute Xi _|_. To this purpose let us rewrite J-'{uj) as follows 

^M=^H+.F(.;')^^^, FM:= (A7) 
If we substitute ()A7|) in the expression for Xi _|_ since = 0, Lemma fA. II implies that 

Hoo , /> ^i([A}—u) )i+i(K —k)u 



,2 , ' C e 

_^ = (cj + \)T{uj ) lim / rfw / c/u;- 



uo 

I III I 



Notice that by introducing a convergence generating factor e~'^'^ with a > in the above 
expression, we can apply the Fubini theorem, and compute the integral over u exactly as in 
Lemma [a. II Hence, we get 

Xi + = (cj' +l)J^{J) lim lim r(l + za') / du -, 

' ° ' ^-.±<x><.-.o+ 1 'J ?^+-(a;2 + l) 

\ ^iii 

[ e-^+M(l,2 + za;z+) ^,(,v^, 

l + ia J ^2 + 1 

Qui 

with i+ = [a~i{K —k)\uq but since the second integral in the above expression is a particular 
case of the second integral entering in ()A5j) . we can conclude that 



_ r(l+ia) i2 /- /" , 
Xi + = ; — [io + \)j-[lo ) lim lim / dio 



u;^ — — 



with Kp = K +ia. Since Lemma FlII. ll implies that the contributions of the frequency intervals 
I± to the wave operators asymptotically at infinity can be made arbitrary small for t —>■ ±00, 
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we extend the domain of integration from Qm to a{Hoo) in the above integral, and compute 
it with the method of contour integrals. Let us define 

C i(u;-J)t 

Xjl^ = lim hm / du , (A8) 

Concerning x{^^ , we complexify the integrand as follows 



^i{z-w )t 

F{z) 



(z2 + l)(y22_^2_^^)lW 

Since for t +00 the imaginary part of e*^* decays exponentially in the complex upper half 
plane, we shall close the contour there. Moreover, F{z) has has two simple poles at i, and two 
branch points at img. We make F{z) single-valued by choosing -di e [0, 27i) and ^2 G (— tt, it] 
with 6*1 := Arg(2; — mg), and ■&2 ■= Arg(2; + rrie) such that \/ — ml = +^yuj^ — ml for 
= 0, and z^ — ml = —-^u'^ — ml for ^2 = vr. To understand how to close the contour 
let us rewrite F{z) as follows 

l+ia 



F{z) 



,i{z-^ )t I J ^2 _ ^2 ^ 



z'^ + 1 \ z^ — ml — 



and analyze the roots of the equation z — mj — k = 0. A simple calculation involving the 
definition of Kp gives = cu''^ —a^+2iaK' . Notice that cj'^— cr^ > ml—a'^+e'^+2mee for every 
uj' E fij. Since we are free to choose a such that < cr < me, it follows that cj'^ — o"^ > 0. 
Taking into account that the sign of an depends on the sign of cu' , we conclude that z'^ lays 
in the first quadrant for e{uj') = +1, and in the fourth quadrant for e{uj') = —1. Hence, it 
follows that z"^ — ml — = possesses for e{uj') = +1 two complex roots , let us say Ci^i,>, 
and C(J2,> = — >, in the first, and third quadrant, respectively whereas it has for e{uj') = —1 
a complex roots co'i^< in the fourth quadrant, and another complex root 0^2, < = —uji,< in the 
second quadrant. 

Let us begin with the case e{uj') = +1. Since we want to apply the residue theorem, we 
choose a contour C such that it circumvents the point ooi^y, and F{z) is analytic within, and 
on C except for the simple pole at z = i. This can be done by fixing the contour C as in 
Figure m Let us define 

gi(aj— |a; \)t 

= I duj 



(cj2 + l){e{uj)^uo'^ - ml - /€p,>)i+^"'> 
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FIG. 1: Integration around the contour C. 



with Kp^> = \k'\ + ia, a'^ = Mml/ \ u' — ml, and F = Fi U U Then, we have 



hm / + hm 

^^0+ I R-^+oo 



V 



/+me-At 



hm / dz F>{z) = 27ri Res(F>(z), z = i) (A9) 

R—>+oo J 



with 



^i{z — \uj \)t 



F>{z) 



(Z2 + l){^Z^-ml - /€p,>)^+*"> 



^i(uj—\u) \)t 



and F3 := [— rrie + /x, me — /I]. Let us first analyze the integrals on C'~\ For i = 1 by 
introducing the parameterization z := m + Jle^^'^ , and taking into account that F^{z) is 
bounded on it can be easily verified that 



j dzFy{z) 


< c/i / 




^0 



vrc/i. 



-(1) 



An analogous relation holds for the integral on C~ . We consider now the integral on C)^' . 
To this purpose let us introduce the parameterization z = Re^^ with < i) < i^Q —7, and 
<?< 'do < 71 /2 where {^q denotes the slope of the ray P in Figure HJ Since F-^{z) is 
bounded on C^"*, we have 



7(1) 





j"do-e 


j dzF{z) 






Jo 
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For R — s> +00 the above integral converges to zero according to Lebesgue dominated con- 
vergence theorem, and the integral on C)j can be treated analogously. Regarding the third 
term in ()A9|) we find for < a < nif. 



+me-IJ. 



< I =<-(me-/x), 



and by taking the limit fi ^ 0^ we obtain 



We consider now 



+me— M r+rue 

lim / duo = I duj ilrrioj) . 



-+me 

lim lim / duiiJu). (AlO) 



Since < l/J(^ml-u;^-\n'\y + a^, and 



it follows that |ilo-(c<j)| < 1/a/c<Jq — m^. Hence, by applying the Lebesgue dominated conver- 
gence theorem to (jAlfljl we end up with 

r-.. ^'r'^""''"'" /M=^^3E4±M. (All) 

Since 2mg — uj'"^ — uj"^ > ml — uj'"^ for \uj\ < m^ it results that 



ml — iu"^ + \k 
mi — uj 



/M < j2 ' UJ enj 



is integrable, and the Riemann-Lebesgue lemma implies that ()A11|) is zero. Regarding the 
computation of the residue at z = +i, it can be checked that 



e * e * 



\Res{Fy{z),z = +i)\ = —— y==- < —j-, \k \ for uenj. 



Therefore, it follows that it does not give any contribution for t +00. Finally, since for 
-R — > 00 the contour — F goes over into the Gamma function contour F' ()A9|) simplifies to 



lim lim Sy = lim lim / dz Fy{z). (A-12) 



r' 
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We reduce now the complex integral in ()A12|) to the Hankel contour integral for the reciprocal 
Gamma function. To this purpose let us rewrite Fy{z) as follows 

with 

F>{z) = ■ , = 

Z-^ + 1 \ Z + UJi^y J 

Since Fy(z) is continuous on, and analytic within F', the integral in ()A12|) is convergent, 
and F^{z) admits a convergent expansion around Wi^y, namely 

oo 

-^>(^) = X^c„(z - u;i,>)" I2: - u;i,>| < ?7, 1] > 0. 



n=0 



Thus, 



, n=0 , 

r' r' 



A simple calculation employing the Hankel contour integral of the reciprocal Gamma func- 
tion (see Erdely et al., Vol.1, p. 14) gives 

/oo 

Since the coefficients c„ depend analytically on a, we can perform the limit cr — O"*", and we 
obtain 

^> r(l + 2a'>)c^'2 

Taking into account that cu' G fi/, we have 

/ . 

lim ^> ~ — TT 5 e'">^°§*, t +00. (A13) 

We consider now the case e{uj') = —1. In analogy to S'> we define 



<x-o+ > r(l + ia'>)cj'2 + i \lj'J 



O-(Hca) 



with /tp^< = — + io", = —Mml/yuj' —ml, and 

gi(2+|a; |)t 



+ 1)(^Z2_^2 _ ft:p_<)l+-< 
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Proceeding as we did for S'>, and taking into account that co'i^< lays on the second quadrant, 
we end up with 

<T^o+ r(l + za^j + 1 / 

Putting together ()A13|1 . and ()A14|1 and taking into account that lu' E fi/, we obtain for 
t +00 

hm / du— \ ~ „ I e*- 



0-^0+ 

cr(//oo) 



(cu2 + l)(eMv/cc^2_^2_^^)iW r(l + la) uj'iu'^ + 1) 
from which it results that 

Let us consider the case t —00. We complexify the integrand function in ()A8jl as follows 

^ pt(Z-Ul )t 

F{z) 



(^2 + l){-^z'^-ml - KpY+i'^ 

Since for t —>■ —00 the imaginary part of e*^* decays exponentially for Im 2; < 0, we have to 
choose the contour in the complex lower half-plane. Moreover, we recall that — m\ = 
— ^uS^ — ml for -^i = 27r, and — ml = +^Juj'^ — ml for ■§2 = — tt. In what follows we 
outline the computation in the case e(ci^') = +1 since the case = —1 is similar. In order 
to apply the residue theorem we fix the contour C such that it circumvents the point — co'i,>. 
Moreover, the function 



^i{z-\ui \)t 



is analytic within, and on C except for the simple pole at —i. The residue theorem implies 
that 

j dz F>{z) = -27riRes(F>(z), z = -i). 
c 

A simple computation shows that the residue is dominated by e*. At this point we can 
proceed similarly as we did for ()A9|) with the only difference that now the part of Fy(z) 
which is continuous on, and analytic within that part of the contour indentating the point 
—uji^y, admits for |2; + co'i^>| < r] with r] > a convergent expansion Yl'^=o'^-n{z + Ui^y)"- such 
that Co tends to zero for a ^ 0^. Repeating the same procedure for e(ti;') = —1, we conclude 
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that Xi^+ ~ for t — oo. Finally, Ti^_ can be obtained directly from Ti^+ by means of 
complex conjugation, and of the transformation t —t. I2- can be computed by the same 
method used for Xi + whereas 12,+ can be derived from X2- by complex conjugation, and 
the transformation t —t. □ 



APPENDIX B: THEOREMS FOR THE EVALUATION OF 



(00) 



Let uj G ri///, oj' G r^/, K be defined as in Section |TT1 k = e{u!')y u''^ — ml, uq > 0, and 
a := Mml/n . Moreover, let T G ^(^m), ^t) be given by (I2H), and 

Mml 



5{t) := —a Log t, a := a(c<j ) = e{uj ) 



' 2 9 



Lemma B.l Let I = I^ := [tu, uo'] for uo < uJ , and I = 1^ := [uj' for uj' < lu, Then, 



U — UJ 



< Mm^ log |t| sup {p{uj)}, p{uj) 



\UJ\ 



{uj'^ — m'lY/'^ 



Proof. The result follows directly from the inequality 



gi<5(t) _ gi5(t) 




/ dx — - — 






Juj dx 



< \uj — UJ \ sup 



duj 



together with the estimate 



duj 



< Mml\og\t\p{uj), p{uj) 



\UJ\ 



□ 

Lemma B.2 

A-t = lim 



P± = lim 

t— >itoo 



du u 



"0 



du u 



itia 



dujj^{uj) (^e^^(*) - e*'^^*) j e'^"^'^ = 0, 

dujT{uj) (e'^^'^ - e*^(*)) ^i{u^-J)t±i{>^' +k)u ^ q 



where the subscript ± attached to A, and P corresponds to the ± entering in the exponents 
of the integrands. 
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Proof. We show the result for A± since P± can be computed with the same method. 
Proceeding as in Lemma fA. II we obtain 



= Mj±i°'r(l±2a') hm hm / dcu J^a)iuj)I^'°'' '^e'^''-'^'^' 

t^±OC (T— »0+ J 

Uiii 

^ y hm hm / rfcu J'm(cu)e-^±M(l,2±za';z±)e*("~'")* (Bl) 

l±ia t^ioo 0-^0+ J 



with J-'(^t){'-^) = J-'{uj){e^^^^^ — e*''*-*''). Let us consider the first integral in (jBlj) . We rewrite the 
integrand as follows 



-77 f, ,\~l^tct -1 x{lu-lj )t 



J^(t){.^) UJ - J K-K ^t{i^-^)t 



uj -uj' K - k' cr =F i{K' - k) [a ip i{K - k)]^^°'' 

where we have used the definition of i± given in Lemma lA.ll Since J-'i^t'j(uj) / (uj — lu') is 
continuous at u; = u' , and has rapid decay for |ti;| — > oo, we can apply the Lebesgue 
dominated convergence theorem to obtain 

:= lim / t/cj J-(t)(cu)(/t-/t')^'"'-'e*("-"')*. (B2) 

t — *iboo / 
Qui 

Without loss of generality let us suppose that e(u;') = +1. For = 1/log^ |t| with |t| > 
which ensures that |ct; | — e > mg + /x we introduce the following decomposition 
of namely 

^iii = {—oo, —rrie — /u] U [rrie + /I, leu | — e'(t)] U [|c<j | — ?(t), \uj \ + e'(t)] U | + e'(t), +oo). 
Let us rewrite e**^'^"'^ as follows 

itdu ' ^ ' 



Concerning the interval / := (— cx3, —rrie — Ji] we use flB3jl to integrate ()B2|1 by parts. In the 
limit t ±oo we have no boundary terms since has rapid decay for |ci;| — » oo, and at 
oj = —me — Jl the corresponding boundary term is dominated by t^^. Hence we get 



Computing the derivative in the above integral, we find the following estimate 

(1), ^ Ci +C2l0g|t| 



A^'^l < hm 
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with constants Ci, C2 > 0. Thus, A^^"* = in the hmit t ±00. Regarding the interval 
// := [me + \lu'\ — e'(t)], we use again ()B3|) . The boundary terms vanish for t ^ ±00 since 
the corresponding boundary term at a; = rrie+Jiis dominated by t^^ whereas at a; = \uj'\—'e{t) 
if we define the function 



we find that 



Thus, we obtain 



hm 



K 



<2^|^(^)| hm 

00 t^±oc 



log 


r2 


\t\ 




t\ 





A 



(1) 



t^itoo it J duj 
II 



d 



lim - / du— (a/c^^ _ ^2 _ 



, \ =Fja -1 



^i{u)~\(jj \)t 



By computing the derivative in the above expression, and by applying Lemma FB. II it can be 
checked that 

|Aff|<C3 hm 1^ 

t— »±oo t 

for some constant C3 > 0. Concerning the interval /// := [\uo'\ —e{t), \uj'\ +e{t)] let 
A^y, = lim / duj J^u){uo){^uo^ -ml- \k e'^^-\''' 

t— >±oo J 

III 

By rewriting the integrand as follows 



and making use of Lemma IB.lf we find that 



IaJV^I < C4 lim 



t^itoo log \ t\ 

for some constant C4 > 0. Finally, for IV := [Icj'I +?(t), +0x3) we define 

a5'^= lim I duo Tu){uj)Wuj^ -ml- |/t'|)^^°'-ie'("-l"'l)*. 

Employing the methods used to compute A^^-*, and A^^ it can be shown that A^y = 0. 
Hence, we can conclude that A*^^) = 0. Concerning the second integral entering in ()B1|) 
since we already showed in Lemma lA. II the boundedness of the function M(l, 2 ± ia ] z±), 
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we can apply the Lebesgue dominated convergence theorem to take the hmit o" — 0^ inside 
the integral, and we end up with the following expression 

A(2) := hm [ rfu;J^(t)(cu)M(l,2±ia';2±) e*("-^')*-^± (B4) 

t — >itoo / 
Qui 

with z± = =i=mo(/t' — k). If we now use ()B3|) to integrate ()B4|) by parts we obtain 

A(2) = -^lirn^l J duj-^ (^J-(t)(cu)M(l,2±m';z±) e"^*) e'^^-^'^\ (B5) 



n 



III 



We do not get any boundary term since J-'{uj) decays rapidly for |ci;| oo whereas at 
uj = ±(me + Jl) the corresponding boundary terms are dominated by t^^. Taking into 
account that 13.4.12 in Abramowitz, and Stegun implies that 

d 



—M l, 2±ia;z±) = ± r - 



M(l, 2±ia'; z±) - ^ M(l, 3 ± m'; z±] 
2±ia 



with M(l, 3 ± za; ; Z-i- + a' where the bound has been obtained by means of 13.2.1 

(ibid.), we can compute the derivative in ()B5|) . and we get the following estimate 



|A(')| < lim 



^5 + C6log|t| 



t^±oo \t\ 

for some constants C5, Cq > 0. Hence, A*^^) = 0, and the proof is completed. □ 
Theorem B.3 Let 

f+oo , p 

h= lim / duu^''^ / rfcu^(cu) e'^w 

J no J 

Uiii 

p+oo , p 

t 5-itoo y 



Uiii 

+00 



i— >itoo 



MO 

fllll 

-00 



>itoo 



MO 



lim / duu-'" / rfcu^(cu) e'^w('^'")+*^^*^ 



with 



V?^^ (cj, ti) = (a; — c<j )t ± (k — K,)u, (p^^^ {u, u) = {u — u )t ± {k + k)u. 
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Then 



h 





ijt- 


^ +00 





ift- 


^ — oo 





ift- 


■» +00 




ift- 


^ — CXD 





ift- 


-> +00 




ift- 


-i> — oo 




ift- 


+00 





ift- 


— oo 



Proof. We give the proof for the first result, the others being similar. By adding, and 
subtracting to e**^*^*-* the term e*^*-*^ we obtain 

Ji = ^ hni^ J^^" du u''^ j dLO J^iuj) (e*^(*) - e^^^) e^'^w^"'") 



'no 



t— >±oo 



The result follows by applying Lemma IB.21 and IA.2I to the first term on the r.h.s. of the 
above expression and to the second term, respectively. □ 

Theorem B.4 Let ip'^^-^{uj , u) , and ip'^^^{uj , u) be defined as in Theorem \B.3l Then 

hoo 



6f = ^lim / duu^'''-' I rfcuj^(cu) e'^w^"'"^^*"^^^ = 



no 



&i = hm / duu^''^-' / rfcu^M e^^w = 0. 



Qui 



Proof. We show the result for &f, the proof for &f being similar. By introducing a 
convergence generating factor e~'^" we apply the Fubini theorem to obtain 

du T{uj) e^(^-^ )*+^^(*) / dT r±*" e~'^\ (B6) 

Ji 



n 



III 
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The integral over r can be computed as in Lemma fA. 11 and ()B6|) becomes 

6f = ±^u^''' lim lim [ duj J^(cj)M(l, l±ia- z±) e^('^-^')*+^'5W-2± 
a t^±oo (7^0+ J 

-M±*"'r(Tia') lim lim [ duo z^"^ e'^'^'^'^'+'^^'\ (B7) 

t^±co 0-^0+ J 

Concerning the second term in the above expression we can immediately take the limit 
0" — s> 0"^ inside the integral whereas for the first term in ()B7|) more care is required since the 
real part of the arguments entering in M coincide. As a consequence 13.2.1 (ibid.) can not 
be used to find a a-independent bound for M(l, 1 ± ia; z±). However, Rez± = mqCt > 0, 
and 13.1.4 (ibid.) imply that for C > there exists a p such that 

\M{l,l±ia;z)\ < \T{l±ia')\ e"«" ( 1 + ) (B8) 



for every z G C\K with K := {z G C| |Re2;| < p, \lmz\ < p }■ Notice that M(l, 1 ± ia; z) 
is bounded for every z E K. Without loss of generality we can choose C = 1, p > 1, 
and ()B8|) gives the estimate |M(1, 1 ± ia;z)\ < 2|r(l ± ia')\ e""*^, from which it follows 
that |e~^±M(l, 1 ± ia'; i±) I < 2|r(l±ia')|. Hence, by applying the Lebesgue dominated 
convergence theorem ()B7j) simplifies to 



ef = ± lim / duj J^(cj)M(l, 1 ± ia; z±) e*(— ')*+i5W--± 

a t^±oo J 



[^mo)^*"r(Tm) lim / dcu J^{iu){K - k)^"^ e 



with z± = =F^(k' — k)uo. Let us define 

6^ := lim / duj J^{uj)M{l, 1 ± ia ; z±) e^(^-^')*+^^W-^±. 

fllll 

By adding, and subtracting a term e**^*^*-* to e**^*^*-* we obtain 
ef = lim [ dtu ^(i)(cj)M(l, 1 ± ia';z±) e*(^-^')*-^± 



►itoo 

Qui 



+ lim e*'^^*) / rfcj J^(cj)M(l,l±ia';z±) e*("^"')*-"± (B9) 

t^±oo J 



II 
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with J-'(_f^{uj) defined as in Lemma lB.21 Concerning the first term on the r.h.s. of ()B9|) we 
use ()B3|) to integrate by parts, and we end up with 

-z± \ i{u}-J)t 



t^±oo it J duj\ ^ '^ ^ ^ ^ 



According to 13.4.12 (ibid.) we have 
d 



M(l, \±ici\z±)T jMil, 2±ia; z±) 

1 ±ia 



— M(l, 1 ± ia; z±) = ±(mo)"^ (- 
duj \uj/ 

which is bounded since M(l, 1 ± za'; z±), and M(l, 2 ± za'; z-|-)are bounded. Hence, we can 

compute S/ by the same method used to evaluate ()B5|) in Lemma fB. 21 and we conclude that 

S/ = 0. Concerning the second term on the r.h.s. of ()B9|) . we apply the Riemann-Lebesgue 

lemma. Hence, ©f = 0. Let us define 

t— >±oo J 
fljlj 

After addition, and subtraction of e**^*-*-' to e*'^^*-' we obtain 



t — J'dioo / 

+ lim e*^(*) / duj J^{uj){k' - k)^'"'' e 

t— >±oo J 



n 



III 



The second term on the r.h.s. of the above expression is zero according to the Riemann- 
Lebesgue lemma. Concerning the other term we use ()B3|) to integrate by parts, and we 
get 

6± = - lim - f du^ (j^u.(uj)(k' - k)^'"'') e'^'^-'^'^' 
t^itoo it J duj \ / 

Qui 

which can be treated by means of the same methods used in Lemma IB. 21 to compute the 
first integral entering in ()B1|) . Hence, Gf = 0. □ 

Theorem B.5 Let h{uj,u) = 0{u~'^) with d^Oiu""^) = 0{u~'^), andw G C^(f2///) such that 
\h{uj,u)\ < w{lo)/u^, and \w' {lo)\ <C{1 + c<j^)" for some constants C, n> 0. Then 

Wf = lim / / rfcuJ^H/i(cu,M) e'^w^"'")+''^(*) =0, 

Juo J 

Qui 



= ^lini^ / duu^''' I dcu^H/i(cu,u) e'^w^"'"^+'^(*^ =0 



UO 



with (p^^.^{uj,u), and ip^_i^^{uj,u) as in TheoreniB. 
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Proof. We show 211^, the proof for 21Jf being similar. Since h(u,iu) = 0{u ^) we can 
immediately apply the Fubini theorem, and we obtain 



2IJf = lim / duj J^{u) e^^(*)Ar(cj)e*^"-'^ >\ N{uj) 

t— »±oo J 

fill 



du M^*" h{uo, u) e 



±i(K —k)u 



uo 



(BIO) 

with |A^(c<j)| < w{uj)/uo. By adding, and subtracting the term e**^*^*^ to e*''^*-', and applying 
the Riemann-Lebesgue lemma, ()B10|) reduces to 

Wf = lim / du J^u^{ij)N{Lj) e'^"^"^'^'. 



Finally, 2nf = follows by using the same method employed in Lemma lB.21 to treat ()B2I 
on the interval I. □ 



Theorem B.6 Let 



Z{uj,u) := < ti-i > ri^[u) 
h{uj, u) 

with h{uj,u) as in Theorem \B.,% and rjoo{u) as in Section 0711 Then 



"0 



iif = lim / duu^''' r]^{u) 



£^ = lim / duu^''' r]^{u) / rfcu J^HZ(u;, m) e'^w^"'"^+'^^*^ = 0. 

t^±oo I I 
J ua J 



Proof. We show the proof for £^ being similar. By applying the Fubini theorem we 
obtain 



^i6{t)(^( \ i{u}-u )t 



£f = lim / du J^{uj) e''^'>Qiuj)e' 



with 



uo 



Q{uj) = / duu riao{u)Z{uj,u) e 

'uo 



±i(«; —k)u 

(Bll) 



\Q{uj)\ < {uq - Uq) sup {\Z\} = {uq - Uq) 



M6|no,woJ 



w{uj)/ul 
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By adding, and subtracting the term e*'''-*'' to 6**^^*^ and applying the Riemann-Lebesgue 
lemma, (jBllj) reduces to 



Sl^ =^lim I did J^(^t){uj)Q{uj) e 



i(u)—uj )t 



Finally, £^ = follows by using the same method adopted in Lemma IB. 21 to treat ()B2|1 on 
the interval I. □ 



Remark B.7 Let us consider the definition of the modified wave operator (jlHj) with 1^ 
replaced by Xqo — X^, i.e. 



s- hm e-^^*(Joo-2:o)e^^-VS~^ 



i— »itoo 



(B12) 



where is such that Too — ioo defines a new identifying operator with cut-off function x 
having compact support := [moj^o], and Uq > Uq > uq > 0. Proceeding as in Section ITTTl 
we can reduce ()B12|) to the computation of the following integrals 

^lini^y" duu^'^^xiu) j rfcu^(cu)Z(M,cu) e'^w("'")+'^^*\ 
lim f duu^"^x{u) I duj J^{u)Z{u,uj) e''^tt)^'^'''^+'^^'\ 

t— >±oo J J 

X njii 

By applying the same method used to prove Theorem IB .61 it can be showed that the above 



expressions are zero. Hence, (jB12|l is zero, implying that our definition (fl3j) does not depend 
on the particular choice of Too- 



Theorem B.8 Let J G ^2//, 1^ = \ ^1 ^ ^' , ot = Mml/ \ ml — u' , and 



1 



Z{uj, u) := < 



,-1 



u 

h{u, u) 



with h{uj,u) as in Theorem W . ^ Then 



lim 

t— >±oo 



"0 

So 



du e-^V j dujT{uj)Z{uj,u) e^(— )*±-«+*'5W = q. 



(B13) 



3^ = lim / dur]^{u)e-'^''u'^ [ du J^{uj)Z{uj,u) 6'^'^-'^'^'^'^''+'^^''^ = 0, (B14) 

t—*±oo I,, J 

fllll 



UQ 
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Proof. By applying the Fubini theorem we obtain 



3f = hm / diu J^{u) e''^'>Q{uj)e' 



+00 



UQ 



with 



IQHI < 



UO 









f \ 


du e"^" < 




> = < 


r(5,xo)r" 











where Xq = f3uQ > 0, and r(-, ■) is the incomplete Gamma function (see Erdely et al., 
9.1.2 Vol.11, p. 136). Notice that for 5 = 1 the incomplete Gamma function gives rise to an 
exponential integral -E'i(xo) = r(0,Xo) which is well defined since Xq > 0. By adding, and 
subtracting the term e*^*-*^ to e*''^*-', and applying the Riemann-Lebesgue lemma, we get 

3f= lim / dLoT^t){uj)Q{uj) e'^^-^'^K 

t — ^iboo / 

flB13|) follows by means of the same method used in Lemma FB. 21 to treat ()B2|) whereas ()B14|) 
can be obtained by proceeding as in Theorem IB. 61 □ 
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